Accidental Degeneracy and Berry Phase of Resonant States 



A. Mondragon and E. Hernandez* 
Institute de Fisica, UNAM 
Apartado Postal 20-364, 01000 Mexico D.F. MEXICO 

a^ 

Abstract 

> ■ 
O 

I We study the complex geometric phase acquired by the resonant states of an open quantum system 

, which evolves irreversibly in a slowly time dependent environment. In analogy with the case of bound 

states, the Berry phase factors of resonant states are holonomy group elements of a complex line bundle 
with structure group C*. In sharp contrast with bound states, accidental degeneracies of resonances 
produce a continuous closed line of singularities formally equivalent to a continuous distribution of "mag- 
I netic" charge on a "diabolical" circle, in consequence, we find different classes of topologically incquivalcnt 

, non-trivial closed paths in parameter space. 

1 Introduction. 
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For many years now, it has been appreciated that there are distinct advantages in describing quantum 
resonances and the quantum phenomena associated to the production, evolution and decay of resonances in 
terms of resonant or Gamow states, since many physical effects are then readily expressed and evaluated Q. 
In this work, we will give closed analytical expressions for the complex Berry phase of an open quantum 
system in a resonant state of a Hermitian Hamiltonian with non-self- adjoint boundary conditions, and we 
will discuss some of its properties. 

During the last fourteen years the geometric phase factors arising in the adiabatic evolution of quantum 
systemsQ have been the subject of many investigations!^, The early literature was mostly concerned 
with the geometric phase factors of closed systems driven by Hermitian Hamiltonians[^. More recently 
there has been a substantial interest in the complex geometric phase acquired by the eigenstates of open 
quantum systems. This problem arises naturally in connection with various experiments which, by their 
very essence, require the observation of the geometric phase in metastable states. The Berry phase in the 
optical supermode propagation in a free laser, which is a classical system described by a Schrodinger-like 
equation with a non-Hermitian Hamiltonian, was studied by Dattoli et al.B. The measurement of the 
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geometric phase in atomic systems with two energy levels, one of which at least is metastable, was also 
described in terms of a non-Hermitian Hamiltonian by Miniatura et al.[Q. The validity of the adiabatic 
approximation for dissipative, two level systems driven by non-Hermitian Hamiltonians was examined by 
Nenciu and Rasche||^], and by Kvitsinsky and Putterman[P], who also established that the Berry phase is 
complex in this case. Sun^] proposed a higher-order adiabatic approximation for two-level non-Hermitian 
Hamiltonians, and showed that the holonomy structure associated to the Berry phase factor in this non- 
Hermitian case is the non-unitary generalization of the holonomy structure of the Hermitian case. 

In two previous papers] 10, 11], we gave explicit expressions for the geometric phase of true resonant 
states, defined as complex energy eigenstates of a Hermitian Hamiltonian which satisfy purely outgoing wave 
boundary conditions at infinity fl^ , and pointed out some of the mathematically interesting and physically 
relevant properties resulting from the extended nature of the singularities in parameter space associated 
with the occurrence of accidental degeneracies of two resonances. In this conection, in another paper we 
showed that the codimension of the accidental degeneracy of n resonances differ significantly from those 
of bound states, we also showed that, close to a crossing of two resonant states, the topological structure 
of the energy hypercomplex surfaces differ significantly from the double conical point singularity typical 
of bound states|^]. Later, by means of a numerical analysis of the experimental data on the 2'^ doublet 
of resonances with T = 0, 1 in the energy spectrum of *Be, we showed, in a realistic example, that a true 
crossing of resonances mixed by a Hermitian interaction may be brought about by the variation of only 
two real linearly independent parameters In this paper we show that the geometric interpretation of 



the Berry phase factor, first given by B. Simon[15| for the adiabatic evolution of closed quantum systems, 
may be generalized, in a very natural way, to the case of resonant states of open quantum systems. That 
is, the adiabatic evolution of resonant states may be interpreted as parallel transport in a complex line 
bundle defined over the space of parameters with structure group C* (the multiplicative group of the non- 
zero complex numbers). Then, the Berry phase factors of resonant states arise as the holonomy group 
elements due to a connection in the bundle such that during the adiabatic evolution the resonant state is 
parallel transported along the fibre. The horizontal spaces are perpendicular to the fibre with respect to a 
generalized inner product of resonant states defined in a rigged Hilbert space. 



2 Resonant States 

2.1 Resonant states in a slowly time evolving environment 

Let us consider the time evolution of a quantum system in a state which is a superposition of unstable 
eigenstates moving in some strong external field of force which changes slowly with time. In order to have 
some concrete example in mind, although a very hypothetical one, we may think of an ^Be nucleus which 
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has only unstable energy eigenstates moving in the field of forces of a double magic nucleus, like ^°Pb, 
in a peripheral collision in which the distance between the two nuclei is never smaller than the sum of 
the nuclear radii. In a semiclassical treatment of the collision, when the centers of the nuclei move along 
classical trajectories, the parameters in the nucleus-nucleus interaction change with time[^]. 

The evolution of the system under the influence of the external perturbation is governed by the time 
dependent Schrodinger equation 

w 

The Hamiltonian H is the sum of the time-independent Hamiltonian Hq describing the evolution of the 
unperturbed system plus a perturbation term Hi which is a function of a number of "external" parameters 
{Xi, X2-..Xj\f} which may change with time, 

H{t) = H^ + Hi{X,{t)). (2) 
The energy eigenfunctions of the unperturbed Hamiltonian are the solutions of the equation. 



Ho^pm {ii , r}j , r) = Sm^m {ii , rij , r) , (3) 

the wave functions ^m{^i-,Vji'^) satisfy the boundary conditions appropriate to a decaying state. 

Assuming that, in the absence of perturbation, the unstable system decays spontaneously in two stable 



fragments, the unperturbed energy eigenfunctions may be written as cluster model wave functions! 17 1 



where (/>^(^i) and (psiVj) the wave functions of the clusters A and B, Umiif) is the radial part of 
the wave function of the relative motion of the two clusters, yjig{r) is a spherical harmonic and A is 
the antisymmetrizer. In our example, ipm would be the eigenfunction of a state of ^Be which decays 
spontaneously in two ^He clusters. In this case Umiir) is a Gamow function. 

2.2 A few facts about Gamow functions 

Gamow functions are the eigenfunctions of the time independent Schrodinger equation which vanish at the 
origen, and behave as purely outgoing waves for large values of the relative distance r, 



Uml{^) = 0, 



(5) 



and 
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^lirn^ [umi{r) - Oi{km,r)] = 0, (6) 
where the function Oi{km,r) is an outgoing spherical wave of complex wave number km and angular 



momentum /[12]. 

The boundary condition is not self-adjoint, in consequence, the energy eigenvalues are complex, 
with Re£m > and Im£m < 0. Hence, Gamow functions are a generalization of bound state eigenfunctions 
in that they belong to complex wave numbers, km = Km — i^m, with Km > 7m > 0, instead of purely 
imaginary ones. This generalization leads out of the Hilbert space based quantum mechanics. Therefore, 
the quantum mechanical rules for normalization, orthogonality and completeness in their usual form, do 
not apply. Nevertheless, bound and resonant states form a bi-orthonormal set with their adjoints, which 
may be extended by a continuum of suitably chosen scattering states of complex wave number to form a 



complete set in terms of which any square integrable function may be expanded |12|. 

The symmetry properties of the Schrodinger equation and the boundary conditions under the operations 
of complex conjugation and reversal of time suggest a definition for the adjoint Umi{km^f) of the Gamow 
function Umi (fcm,r)p^, 

Uml{km',f)=Uml{ — km)f)- (J) 

Now, if the adjoint Umiikm^f) is identified with the bra-eigenfunction, the quantum mechanical inner 
product, or bra-c-ket rule, may be generalized to include Gamow eigenfunctions. With this prescription, 
matrix elements of operators and inner products are computed in configuration space representation as 
integrals over the radial variable r. Since the Gamow function Umiikmif) and its adjoint Umi{km, r) oscillate 
between envelopes that grow exponentially with r, the integrals over r must be properly defined. This may 
be done, either by analytic continuation in the complex /c-plane from above ^ or by means of an Abel 
regulator|l£] with a suitable convergence factor and a limiting procedurepO, 21 1. Both procedures give the 



same result|22, 23|. In this paper, we will adopt the second definition. 



Then, it may be shown that 
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lini / e ^''^ u*^i{-k*^,r)unv{kn,r)dr = Smn^iu (8) 



i.e, the Gamow eigenfunctions are orthonormalized in a generalized sense. It may also be shown [12, 24 1 
that, for any two square integrable functions, /(r) and g{r)^ the following relation holds 

S^nr)g{r)dr = Y.n{U^nr)unl{kn,r)dr] x [^^ uU-kiy)g{r')dr']} + 
jjk [J,^ r{r)Mk,r)dr] [J,^ {-k* ,r')g{r')dr'] . 
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In this expression, the functions u^i^r) are bound states or Gamow state eigenfuctions belonging to real 
negative or complex eigenvalues the functions (j)i{k,r) are scattering partial wave functions of angular 
momentum / and complex wave number k. The integration contour C, in the wave number plane k, starts 
from the origin as a straight line with slope -1, it goes down to Imk = a and then, it continues as a 



straight line parallel to the real axis 1 12]. Finally, the square brackets around the integrals mean that, when 



necessary, the integrals are defined by means of a gaussian regulator and a limit as in (M) or by analytical 



continuation from above as in Mondragon and Hernandez 1 12 



Since /(r) and g{r) are arbitrary functions, we are justified in writing the expansion 

5(0 = ^Umiir) < Umi\g > + J^(j)i{k,r) < (l)i{k)\g > dk, (10) 



m 

the index m runs over bound and resonant states. 
The expansion coefficients are given by 



^2 



< Umi\g >= lim / e u*{-k*^,r)g{r)dr, (11) 

< Uk)\g >= lim r e~'''"<PU-k\r)gir)dr. (12) 

Once the validity of the quantum mechanical inner product has been generalized to apply to bound 
and resonant eigenfunctions of the relative motion of the two clusters, the generalization of (P) and (^) to 
expansions of many body wave functions in terms of cluster model bound and resonant eigenfunctions is 



fairly straightforward] 12, 24 1. 



2.3 The mixing matrix 

We may, now, go back to our problem, namely the time evolution of the resonant states of a many body 
system moving in a slowly time-dependent external field of force. Since we are interested in the time 
evolution of a state ^ which is a superposition of unstable states, we make an expansion of the wave 
function ^ in terms of bound and resonant states of Ho, 

= J2amit)ip^{^i,Vj,r)+ f b{k;t)ip^+\k;Ci,rjj,r)dk. (13) 

In general, the index m runs over bound and resonant states. In our example, (fm would be the complex 
energy eigenfunctions of the ^Be nucleus which decays spontaneously in two ^He clusters. The scattering 
states f^^\k, Ci,r]j,r) of complex wave number k and the integration contour C are defined in the previous 
subsection. 

Substitution of (IS) in (|l|) gives the set of coupled equations 
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dt 



-^£ma„,{t) -^J2< y^rn\Hi{t)\ipn > an{t) " ^ / < y^rn\Hi{t)\ip'^+^k) > b{k;t)dk (14) 



and a similar expression for db{k,t)/dt. 

We will use the notation IfmidiVji''') > fo^' the Gamow function, and < ifmiS,i,r]jj''")\ for its adjoint. 
Hence, the matrix element of the perturbation term Hi{t) taken between bound or resonant states of the 
unperturbed system is given by 



\Hi{t)\ ifn) = J ■■■ J {¥'m{^i,Vj,r)\Hi{t)\ipni^i,Vj,r))d^^i ■ ■ ■ d-^U-id-^Vi ■ ■ ■ d-^ijB^id^r, (15) 

the integral over the radial variable r is defined as in (|^. 

When the interactions are time reversal invariant, the dual of the complex Gamow function Umiir) is the 
same function 1 12, 24 1. But, when the interactions are not time reversal invariant, the Gamow function 



and its dual are not the same function. 

Therefore, when the forces acting on the system are time reversal invariant, the complex matrix H, with 
matrix elements 

Hmn(i) = £mSmn + {fm \Hl{t)\ (fn) (16) 

is symmetric but non-Hermitian. When the forces acting on the system are not time reversal invariant, H 
is, in general, complex, non-symmetric and non-Hermitian. 

The contribution of the non-resonant background integral over the continuum of scattering wave func- 
tions to ( p!3[ ) and (^) will not be relevant to the following discussion. Therefore, to ease the notation, we 
will disregard the contribution of the background of scattering functions. With this truncation. 



m 

Then, the set of coupled equations (|l^) reduces to 



(17) 



.(t) = -^EHmn«n(i), (18) 



where am{t) is the time derivative of am{t). 

3 Geometric Phase of a Resonant State. 



The time evolution of the quantum system is governed by the Hamiltonian or the matrix H, which are 
functions of N real, linearly independent parameters. Therefore, we may consider the matrix H embedded 
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in a population of H matrices smoothly parametrized by N external parameters which take values in some 
domain D of a manifold or parameter space. Each point in D represents an H matrix. When the external 
parameters change with time, the system traces a path C in parameter space. In the following, we will 
study the behaviour of the system in the time interval [0,T], assuming that at the initial time, t = 0, it is 
in an eigenstate of H, at t = T the parameters {Xi,X2, ....Xn) have returned to their initial values, and 
the adiabatic theorem holds. Then, the unstable system traces a closed path C in parameter space while it 
remains in an eigenstate of H{t) at all times. 

In the absence of symmetry, H has no repeated eigenvalues at almost all points in the domain D. The 
set of points in D where H has one twofold repeated eigenvalue is a subdomain D' C D. If D has N 
dimensions, and D' has N' dimensions, then N = N' + k, k is the codimension of the twofold degeneracy. 
Similar relations hold in the case of an m-fold degeneracy. 

It will be assumed that the complex non-Hermitian matrix H has no repeated eigenvalues at all points 
on the path C. In consequence, at all points on C, it may be brought to diagonal form by means of a 
similarity transformation 

K-^HK = E, (19) 

where E is the diagonal matrix of the complex energy eigenvalues. The columns in the matrix K are the 
instantaneous right eigenvectors of H. In an obvious notation 

K=(|<^«),|<^(2)^^ |<^W)....|0H)). (20) 

which satisfy the eigenvalue equation 

H{Xi{t))\^^^\t))=m\<P^'Ht)). (21) 
The rows in are the corresponding left eigenvectors of H, properly normalized, 

{^^(i)\cl>U)) = 6ij. (22) 

The adiabatic basis, {\ips {£,i,r]j,r; Xi{t)))} , of instantaneous bound and resonant energy eigenstates of 
the complete Hamiltonian H, is obtained from the set of unperturbed bound and resonant states with the 
help of the matrix K 

\'Ps{^i,Vj,r;Xi{t))) = '^\(pm{^i,r)j,r))Kjns{t) (23) 

m 

and their adjoints (duals) are given by 
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{^s{^„Vj,r;Xim = J2{K-\t))sn{M^i,Vj,r)\. (24) 

n 

It follows that the instantaneous bound and resonant energy eigenstates of H(t) satisfy an orthogonality 
relation similar to (^), which they inherit from the unperturbed complex energy eigenfunctions 

{'fs{^i,rij,r;Xi{t))\(fs'{^i,Vj,r;Xi{t))) = J2m,n{'^'^)sn{^n{^i,Vj,r)\(pm{Ci,Vj,r))Kms' = 

_i 

If we multiply both sides of eq.(|23|) by K~^, we obtain 

\y^rn{^i,r]j,r)) =^\ips{^i,Vj,r;Xi{t))){K~'^)sm- (26) 

s 

Substitution of this expresion in (^), gives the expansion of ^' in instantaneous energy eigenfunctions 
of H, 



^ = E l'^- iC^,Vj,r■,Xl{t)))a,it), (27) 

s 

where 

a,(t) = V(K-i(t)) a„(t). (28) 

^ — ' \ / sn 

n 

In (|27| ) we have kept only the summation over bound and resonant states and, as in the previous section, we 
have disregarded the contribution from the integral over the continuum of scattering functions of complex 
wave number. 

Similarly, the expansion of becomes 

H^! = Y,\^s Hi ,Vj,r;Xi{t)))t^s {t)a, (t). (29) 



Substitution of (^) and ( p9[ ) in the time dependent Schrodinger equation gives the set of coupled 
equations 

+ V {ipslVsfiipra) • ^am{t) = -i^s{t)as{t) (30) 

at ^ at 

m=l 

where, |Vs({(^m > is the gradient of \ipm > in parameter space, 

= E [gjrj^m >j ^ = W^^m > (31) 

It will be assumed that the non-adiabatic transition amplitudes are very small 
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1 

la.. 



<< 1, m^s. (32) 



Then, we can make the approximation 



^-^-^m-{^s\V^^s)■^. (33) 



Integrating both sides, we get 



as{t) = exp[-^ f ^s{t')dt']exp[i^s]as{Q), (34) 
ri J to 

the first factor is the complex dynamical phase, whereas the second one is the complex Berry phase given 
by 

ls = i J ((^slVsRc^s) • d§:. (35) 

Direct evaluation of 'Vsalifs) requires a locally single valued basis for \{ps) and can be awkward. Such 
difficulties are avoided by transforming the path integral into a surface integral with the help of Stokes 
theorem 

ls = iy2 / (VsR<^5|<^m) A ((^m|VsR<^s) • dT,, (36) 

where E is any surface in parameter space whose boundary is the curve C, and A means wedge product. 
This expression may be written in a more convenient form by means of the identity 

{(Psl^^i'fm) = A ^_ A {'fs\^^Hi\(p,n), (37) 

then, 

^s = i 2^ / -T — -rT2 ' (3^) 

^^7e79E=c (4 -An) 

provided the surface does not cross any point in D where the denominator vanishes. Since the dependence 
on |Vs({99s) has been eliminated in (^), phase relations between eigenstates with different parameters are 
no longer important, any complete set of instantaneous eigenfunctions of the time dependent Hamiltonian 
H may be used to evaluate the integral in (|38[). 

Explicit expressions for the Berry phase in terms of our parametrization of the interaction Hamilto- 
nian may easily be obtained. Since the unperturbed Hamiltonian and its bound and resonant states are 
independent of time, the time dependence of the bound and resonant instantaneous energy eigenstates of 



9 



the perturbed system is entirely contained in the matrix K, which is a function of t through the time 
dependence of the external parameters. Therefore, from eqs.([2^), ( p^ ) and ( p5[ ) we get 



Is = « 



^ [k-^(VjjK) 



(39) 



The columns in VjjK are the gradients of the instantaneous eigenvectors of H(t). Therefore, 7^ may 
also be written as 



(40) 



Furthermore, from (SC) we may derive a sum rule for the geometric phases of the interfering resonant 



states. Writing the integrand in the right hand side of (3£) in terms of K and VsrK, and taking the sum 
over s, we get 



^ls=i / tr 



(V^K-i) A (VjjK) 



(41) 



but 



1 

— tr 

271 



(VjjK-i) A (VjjK)] =Ci{N), 



(42) 



is the first Chern class of a complex vector bundle defined by the matrix H over the parameter space M. 
Its integral is the first Chern number [B^] 



1 

2tt Jt. JdT, 



tr 



(VsrK-1) a (VsrK) 



• dS = ci. 



(43) 



Hence, 



Y,ls = iftr \k-^ (VsrK)] • d^. = -27rci 



(44) 



is a topological invariant. 



4 Berry phase factors of resonant states and holonomy in a complex 
line bundle 

4.1 Resonant states as elements of a rigged Hilbert space 

Since 1983, B. Simon p5[| pointed out that for Hermitian closed quantum systems, the adiabatic evolution 
can be interpreted as a parallel translation in a Hermitian line bundle and the Berry phase factor is the 
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holonomy in such a bundle. In this section , it will be shown that this geometric interpretation may be 
generalized in a very natural way to the Berry phase factor of resonant states of open quantum systems. 

The resonant or Gamow functions {\4>mi(,i, Vjj''" >} of the unperturbed system are eigenfunctions of the 
self-adjoint unperturbed Hamiltonian Hq, which satisfy purely outgoing wave boundary conditions for large 
values of the separation distance r of the decay fragments. Since the boundary condition is not self-adjoint, 
the corresponding energy eigenvalues Sm, are complex, with ReSm >0 and Im£m <0. In configuration 
space representation, the Gamow functions, as functions of r, behave as outgoing waves which oscillate 
between envelopes that increase exponentially. In consequence, Gamow functions are not square integrable 
and cannot be characterized as elements of a Hilbert space. 

Gamow functions are usually associated with the resonance poles of the scattering matrix and the re- 
solvent operator of the time independent Schrodinger equation which lie in the lower half-plane of the 
unphysical sheet of the Riemann energy surface In order to give a proper mathematical charac- 

terization of Gamow functions as elements of a space, one has to specify a rigged Hilbert space in which 
the Gamow states {|93m(Ci) ''') >} are defined as generalized eigenvectors of the Hamiltonian Hq with 
generalized complex eigenvalues £m- 

Following Bohm and Gadella||2^], we will associate to the resonant poles of the resolvent operator of the 
Schrodinger equation of the unperturbed system, a rigged Hilbert space 



in this expression, is the space of well-behaved functions of the position coordinates {(,i,rjj,r) which 
are Hardy class 2 functions of the complex energy E from above, 7i is the Hilbert space of square integrable 
functions and <I>^ is the space of continuous antilinear functionals on In this way, the Gamow states 
of the unperturbed system are defined as continuous antilinear functionals on <I>-|-, that is. 



(45) 




(46) 



The adiabatically evolving Gamow functions \(ps{S,i,rjj,r;^{t)) > of the complete Hamiltonian H, in- 
troduced in section 3, eq. ([2^) , are linear combinations of the Gamow functions l^mi^iTVj^''') > of Hq, 




(47) 



m 



where K[5?(t)] is the matrix which diagonalizes the complex non-Hermitian matrix H[3f?(t)]. 



Therefore, the adiabatically evolving Gamow functions are also elements of <I>1 



(^,(e„r/,-,r;[5}?(t)]) >e$; 



(48) 
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The dual of the Gamow eigenvector |(^s[JR] > is the eigenvector < = J2m(^~^)sm < ^m{CiiVj)\ 

defined in eg. (p^) , corresponding to the same complex eigenvalue <?s[^]- Since, by assumption H[5R] has no 
repeated eigenvalues for 3?eC C M, the adiabatically evolving Gamow eigenvectors and their duals satisfy 
the orthogonality relation 

< ifsm^s'm >= 0, Ss^s's (49) 

inherited from the Gamow eigenfuctions \^Pm{£,i,Vj:''^) > of Hq, and they may be normalized to one 

<^smv>sm>=i (50) 

Having characterized the adiabatically evolving Gamow functions as elements of a rigged Hilbert space 
we may turn to the question of the geometric interpretation of the Berry phase factors of resonant states. 

4.2 Adiabatic evolution and parallel translation 

Let us suppose that, as time varies, the self-adjoint Hamiltonian H[?R.{t)] and its instantaneous energy 
eigenstates | r^j, r; 5}?(t)) > make a cyclic excursion in a closed circuit C in parameter space. We 

assumed that the time dependence of the Hamiltonian justifies the adiabatic approximation. We will 
further assume that, if 3^(t) is any point on C, the complex matrix H[3f?(t)] has no repeated eigenvalues 
in an open neighbourhood of 5Re C, that is, the Hamiltonian H(t) has no repeated instantaneous complex 
eigenvalues when 3f?e C, and that both the generalized eigenvectors | (^^(^j, r/j, r; 5R(t)) > and the eigenvalues 
iSs[5R(t)] are smooth functions of 5ReC. Moreover, we will also assume that the evolution of each state is such 
that there are no level crossings along C. 

To avoid a clumsy notation in the rest of this section we will not write the position coordinates {(,i,rjj,r). 

Let us call Fg the complex line bundle defined by the Gamow eigenfunction | (^s[5R(t)] > over the 
parameter space M of the adiabatically evolving system 

Fs = {{^,\ipsm >)\H[^]\ipsm >=£sm\0sm >,§eM} (si) 

Its fibre is a complex, one-dimensional linear space 

L|j := [li^sm > >=e'''\^sm >,\^sm > e<!>l,ameC*MeM} (52) 

where <1>^ is the space of antilinear functionals defined over the space of well-behaved functions of 
Hardy class 2 from above and a [5ft] is a complex function of SfteM. 
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Under the assumptions made above, the quantum number s, labehing the eigenvalue and the 

eigenstate \tPsW\ > of the instantaneous Hamiltonian, is an adiabatic invariant. Hence, a system prepared 
in a state |^s[5R(i)] > such that 

\i>smto)] >= Csmomsmto)] > (53) 

will evolve with H[^{t)] and be in a state 

mm] >= CsimUsim] > m 

at t. 

Now, consider the decomposition of a tangent vector in vertical and horizontal parts 

|l4W>=(|l4w)^^ + (|4(t)>)^ (55) 
The horizontal part, orthogonal to the fibre is 

J^lMt) >) = E < 'Pmrnij^lMt) > \0mm > (56) 

and the vertical part, along the fibre is 

>)|^ =< ^smij^iMt) > \^sm > (57) 

If IV's(t) > is an evolution state in an adiabatic change, 

j^lMt) >) ^1 = < hatcps[^]\j^0s[m] > Cs[^]j I'fsi^] > (58) 

The condition for parallel translation along the curve C is 

<(^,[5R]|^4[5ft] >=0, SReC (59) 

which gives a one-form equation 

dCs m +<^s m I d0s [^] > Cs =0 (60) 
Its solution gives the Berry phase factor of the resonant state 

exp[i I i < (ps\d hatips >] = exp[z7s[5R]]. (61) 

J c 

In the case of a cyclic evolution, such that 3ft(0) = 3?(r) and C is a closed curve, the complex phase 7^ 
may be expressed as 

13 



-fs{C)=i£<ips\dips> ■ (62) 

Now, it will be shown that the condition for the adiabatic evolution of a resonant state is equivalent to 
the condition for parallet translation along C. 

The condition for the adiabatic time evolution of a resonant state, equation (^), may be written as 

e^/o {I (a.(t)e^/o'^»(*')'^*') + < ^.||^. > a.(t)e^ /o = 0. (63) 

Now, making use of the normalization condition, eq. (|25|) , this expressions may be written as 

e-^r ^^(*')-' < [mm)] > ^ + ^^*^C.(t)} = 0, (64) 

where Cs{t) is given by 

C,(t) = e^/o^^(*')'^*'a,(t), (65) 

multiplying both sides of (U) by the dynamic phase factor, the condition for adiabatic evolution of a 
resonant state takes the form 

< ipslj^i^s >= 0, §eC (66) 

which is the condition for parallel translation of the resonant state |V's[3^(^)] >= C's[3^(0]l'^s[3'?(i)] > along 
the curve C. 

It follows that, except for the effect of the dynamical phase factor, the condition of adiabatic evolution 
of a resonant state is equivalent to the condition of parallel translation on the complex line bundle Fg. 

The geometric phase factors exp[i^sic)], occurring in the adiabatic evolution, are holonomy group ele- 
ments of the complex line bundle Fg . 

5 Accidental degeneracy of resonances. 

As is apparent from (|38| ) , non-trivial phase factors of the energy eigenvectors or eigenfunctions are related to 
the occurrence of accidental degeneracies of the corresponding eigenvalues. In the absence of symmetry, de- 
generacies are called accidental for lack of an obvious reason to explain why two energy eigenvalues and £m 
of H should coincide. However, if the matrix H is embedded in a population of complex non-Hermitian ma- 
trices [il{Xi, X2, ...X]\f)] smoothly parametrized by N external parameters {Xi,X2, ....Xn), degeneracy in 
the absence of symmetry is a geometric property of the hypersurfaces representing the real or complex eigen- 
values of H in a (A^+2)-dimensional Euclidean space with Cartesian coordinates {Xi,X2, ....Xn, Re£, Im£). 



14 



The energy denominators in (|3^ ) show that when the circuit C hes close to a subdomain D' d D, va. pa- 
rameter space, where the matrix H has an m-fold repeated eigenvalue, and the state \(ps) is involved in this 
degeneracy, the Berry phase 7s (c) is determined by the geometry of the energy hypersurfaces close to the 
crossing of eigenvalues and the other (m-1) states involved in the degeneracy. 

It is in this connection of the Berry phase with the accidental degeneracy of energy eigenvalues that the 
non-Hermiticity of the matrix H plays an important role. In contrast with the case of Hermitian matrices, 
square, complex, non-Hermitian matrices with repeated eigenvalues cannot always be brought to diagonal 
form by a similarity transformation. However, any n-dimensional, square complex matrix H may always 
be brought to a Jordan canonical form E by means of a similarity transformation. 

If H has V [u < n) different eigenvalues, ....E^^ with multiplicities ^i{Ei), the Jordan canonical 

form E is the direct sum of v square Jordan blocks Ej. Each Jordan block Ej is the sum of a diagonal 
matrix if^jl^.x/^j, and a nilpotent matrix N^-. Corresponding to each Jordan block, there is a cycle of 
lenght Hi of generalized eigenvectors. When the lenght /Xj of the cycle is >2, the codimension of the 
accidental degeneracy, the geometry of the energy hypersurfaces at the crossing and the properties of the 
generalized eigenvectors involved in the degeneracy of resonant states differ substantially from those of 
bound states ( []T3|| ). Rather than trying to develop a theory of the most general case, in the following we will 
examine the simplest possible case, namely the accidental degeneracy of two resonances and the topology 
of the energy surfaces close to a crossing of two resonances in parameter space. 

6 Degeneracy of two resonances 

Let us consider a system with two resonant states strongly mixed by a Hermitian interaction, all other 
bound or resonant eigenstates being non-degenerate. We may suppose that we already know the correct 
eigenvectors of H for all the real and complex eigenenergies Sg, except for the two the crossing of which 
we want to investigate. Using for this two states two vectors which are not eigenvectors but which are 
orthogonal to each other and to all other eigenvectors, we obtain a complete basis to represent H. In 
this basis, H will be diagonal except for the elements H12 and H21. The diagonal elements Hn and 
H22 will, in general, be non-vanishing and different from each other. There is no loss of generality in 
this supposition, since any complex matrix H may be brought to a Jordan canonical form by means of 
a similarity transformation. When the eigenvalues are equal, H2X2 is either diagonal or equivalent to a 
Jordan block of rank two. Hence, we need consider only the conditions for degeneracy of the submatrix 

H2x2- 

The matrix H2X2 may be written in terms of the Pauli matrix valued vector a = (cri, cr2, cjs) and the 
2x2 unit matrix as 
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H2X2 =£1+(^R- i^f^ ■ a, (67) 

where R and F are real vectors with cartesian components (Xi, X2, X3)and (ri,r2,r3). When the forces 
acting on the system are time reversal invariant, X2 and r2 vanish. 

In the absence of more specific information about the external parameters X, we will parametrize H2X2 
in terms of R and T. Then, 

i: = R-i^f. (68) 



From (67), the eigenvalues of H2X2 are given by 

Ei,2 = ^^Te. (69) 



where 



R-i^f^ . (70) 

Then, Ei and E2 coincide when e vanishes. Since, real and imaginary parts of e should vanish, the 
condition for accidental degeneracy of the two interfering resonances may be written as 

Rl - = 0, (71) 

and 



Rd-fd = 0. (72) 

These equations admit two kinds of solutions corresponding to H2X2 being or not being diagonal at the 
degeneracy: 

i) When both and vanish, eqs. ( [7l| ) and (^) define a point in parameter space, H is diagonal at the 
degeneracy and the submatrix H2X2 has two cycles of eigenvectors of lenght one each. The vanishing of 
implies that Im£ also vanish [Q, therefore the two complex eigenvalues of H which become degenerate fuse 
into one real positive energy eigenvalue embedded in the continuum. Since in this case all the Cartesian 
components of Rd and F^ should vanish, the minimum number of linearly independent, real, external 
parameters that should be varied to produce a degeneracy (codimension) of two resonances to form a 
bound state embedded in the continuum is four or six depending on the quantum system being or not being 
time reversal invariant. 

ii) In the second case, when the degeneracy conditions (Iti] ) and ( ^2|) are satisfied for non-vanishing 
Rd and F^, these equations define a circle in parameter space. In this case H2X2 is not diagonal at the 
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degeneracy, it is equivalent to a Jordan block of rank two and has one cycle of generalized eigenvectors of 
lenght two. In this case, the two complex eigenvalues become one two-fold repeated eigenvalue 8. Since, 
the two linearly independent conditions ( [7l| ) and ( [7^ ) should be satisfied for non-vanishing values of and 
Td, at least two real, linearly independent parameters should be varied to produce a rank two degeneracy 
of resonances. Hence, the codimension of a degeneracy of resonances of second rank is two, independently 
of the time reversal invariance character of the interactions. 

6.1 Computation of the geometric phase. 

Let us consider now the computation of the Berry phase in the simplest situation when the degeneracy 
involves only two resonant states. The matrix 7i that mixes the two interfering resonances off degeneracy 
is 



n 



^ Z -i\T X-iY 
yX + iY -Z + i\T 



(73) 



This matrix has two right and two left eigenvectors and may be diagonalized by a similarity transformation 

/ 



-e 
e 



(74) 



where 



K=^\ I. (75) 



^ and Tj are short hand for X — iY and Z — i^T respectively.The matrices K^^ and V^jK are readily obtained 
from (^0|) and ([75|). If we assume that the interfering resonances are mixed by a Hermitian interaction, T 
will be kept fixed. Then a straightforward calculation gives 



1 r 1 



71 



2 7c re (e + 7?) 
and 



(f X ^) • dR, (76) 



These expressions are very similar to the well known results obtained for the geometric phase of bound 
states An obvious difference is that the geometric phase of resonant states is complex since e and r] are 
complex functions of the parameters R and T. There is another important but less apparent difference: In 
the case of an accidental degeneracy of resonances (F 7^ 0), the denominator in the right hand side of (^) 
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and ( |77|) vanishes on the continuous hne of singularities defined by eqs.(^) and (|7^), which will be called 
the diabolical circle, and not at one isolated point as is the case for bound states. It follows that two kinds 
of non-trivial, topologically inequivalent closed paths are possible. First, those paths which surround the 
diabolical circle but are not linked to it. Second, the closed paths which are linked to the diabolical circle. 
Paths of the first kind are clearly analogous to the non-trivial paths that go around the diabolical point 
of bound state degeneracies while paths of the second kind have no analogue in accidental degeneracies of 
bound states. 



Diabolical 
Circle 




Figure 1: In the evaluation of the Berry phase of two interfering resonant states there are two kinds of non- 
trivial, topologically inequivalent closed paths in parameter space. First, those, like C^^\ which go around 
the diabolical circle but are not linked to it. Second, those, like C^^^\ which turn around the diabolical 
circle and are linked to it. 



For paths of the first kind it is always possible to find a surface E which spans the closed path C and 



does not cross the diabolical circle. Then, Stokes theorem applied to ( [Tq) and (77) gives 



Is 



^\ 2 



(78) 



where s = 1, 2. Since 72 changes into 71 when Ei and S2 are exchanged and the sign of dS is changed, the 
normals for S2 and Si should be oppositely oriented. If we say that F points upwards, then Si is above 
and S2 is below C. 

This is the same result as would have been obtained from the general expression (^8|) and our parametriza- 
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cion of the perturbation term in the Hamiltonian; no summation over intermediate states occurs in (^) 
since in the simple case of only two interfering resonant states, the summation in (^8|) has only one term. 
Adding 71 and 72, we get 



71 +72 



(79) 



where S is a sphere with the diabolical circle completely contained in its interior. The integral is easily 
computed when R > ^T, the result is 



-2^(71+72) = !, (80) 

in agreement with our identification of the integral in ( |4^ ) as the first Chern number of the complex line 
bundle defined by the eigenvector of the non-Hermitian matrix H2x2- 

It is now easy to show that the resonance degeneracy produces a continuous distribution of singularities 
on the diabolical circle. The surface integral in (^) may be written as a volume integral by means of Gauss 
theorem. Then, 



71 +72 



V 



{R-^^T 



dV, 



(81) 



where V is the volume inside J2 ^-nd bounded by it. The term in round brackets under the integration 
sign vanishes when - i^r) / 0. Therefore, the non-vanishing value of 71 + 72 implies the occurrence 
of (5- function singularities of the integrand on those points where ^R — ^^T^ vanishes. 
Hence, 



R-ihT 



R-i^r 



s{R-^r 

R^ 



6 (cos 9) , 



(82) 



the factor multiplying the delta function is needed to reproduce the value 2tt of 71 +72- 

Let us turn to the case of closed paths of the second kind, that is, those paths which are linked to 
the diabolical circle. In this case there is no surface S which spans the closed path c^^ without crossing 
the diabolical circle. Therefore, it is not possible to use Stokes theorem to convert the path integral into 
a surface integral. However, we may still compute the geometric phase from the path integral. To this 
end, we change from Cartesian coordinates {X, Y, Z) with OZ parallel to f , to spherical coordinates in 
parameter space. Then, eqs. (76) and ([77]) take the form 
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I r I r ( R COS 9 — {■^Fj dip 



(83) 



^ ^R'^ -\T^ -iVRcosO 



the path c is specified when R and 9 are given as functions of ip. 

From (p3|), it follows that, for closed paths which are linked to the diabolical circle, 



7i + 72 



/ = 0, 



c- ' of secondkind 



(84) 



since, in this case, the angle ip starts out at some value ipQ and, as the system traces the path c , it 
oscillates between a minimum and a maximum values and finally ends at the same initial value 999 • There 
is no analogue to this case in the geometric phase of bound states. 

7 Results and conclusions 

The purpose of the foregoing has been to discuss the geometric phase acquired by the resonant states when 
they are adiabatically transported in parameter space around a degeneracy of resonances . 

As in the case of bound states, the condition of adiabatic evolution of resonant states may be given a 
geometric interpretation as parallel translation in a complex line bundle Fg. The geometric phase factors 
acquired by the resonant states when adiabatically transported in a closed circuit in parameters space are 
holonomy group elements of the complex line bundle Fg. 

In the case of two resonant states mixed by a Hermitian interaction we find two kinds of accidental 
degeneracies which may be characterized by the number and length of the cycles of instantaneous energy 
eigenfunctions at the degeneracy. In the first case there are two linearly independent eigenfunctions be- 
longing to the same real positive repeated energy eigenvalue, that is, two cycles of lenght one. In the 
second case there is only one resonant eigenstate and one generalized resonant eigenstate belonging to the 
same degenerate (repeated) complex energy eigenvalue, i.e. one cycle of lenght two. At degeneracy, the 
Hamiltonian matrix has one Jordan block of second rank. 

In the generic case of a non-time reversal invariant system, when the degeneracy is of the second rank, 
the topology of the energy surfaces is different from that at a crossing of bound states. The energy surfaces 
of the two resonant states that become degenerate touch each other at all points in a circle. Close to the 
crossing, the energy hypersurface has two pieces lying in orthogonal subspaces in parameter space. The 
surface representing the real part of the energy has the shape of a hyperbolic cone of circular cross section, 
or an open sandglass, with its waist at the diabolical circle. The surface of the imaginary part of the energy 
is a sphere with the equator at the diabolical circle. The two surfaces touch each other at all points on the 
diabolical circle [^]. 
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Figure 2: Two interfering resonances whieh initially have equal level energies but different half widths are 
mixed by the Hermitian interaction R- a. As R increases from to ^F, the points representing £i and £2 
approach each other along a meridian circle on the sphere representing Ime in parameter space, until they 
meet at a point on the equator, which is the diabolical circle. At this point the two resonances become 
degenerate and the mixing matrix is equivalent to a Jordan block of rank two. When R becomes larger 
than ^F the points representing £1 and £2 move away from each other on the hyperbolic cone representing 
Ree. 

In the case of two interfering resonant states, the geometric phase acquired by the resonant states when 
transported around the diabolical circle in a closed path which is not linked to it, may be written as the 
sum of two terms. 

The first term, ^y^'^'^d- ^c^^ ^ is the real geometric phase which a negative energy eigenstate would have 
acquired when transported around a diabolical point in a closed path in the same parameter space. The 
second term is complex, it gives rise to a change of the phase and a dilation of the resonant state eigenfunc- 
tion. Its imaginary part may be positive or negative, in consequence, it may produce an amplification or 
a damping of the wave function which may compensate or reinforce the attenuation due to the imaginary 
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part of the dynamical phase factor. 

When the resonant states are transported in a closed path c^^ which does not go around the diabolical 
circle but is linked to it, the geometric phase they acquire is 



7i7(c^^)=±A7(c^^). (86) 

Since it is not possible to find a continuous surface Yl, which spans the closed path c^^ without crossing 
the diabolical circle, we can not make use of the theorem of Stokes to convert the path integral into a 
surface integral. However, it may readily be computed as a path integral from the expression 

A7 C^0=/ \ . ^ , (87) 



which is obtained from (83) and (p^). As in the previous case, A7 (^c^^^ is complex and produces changes 
of phase and dilations in the resonant state wave function. This case has no analogue in bound states. 

The sum of the geometric phases acquired by two interfering resonant states which are transported 
around a degeneracy in a closed path of the first kind in parameter space is a topological invariant, namely 
the first Chern number Its value is the "magnetic charge" on the diabolical circle. For paths of the 
second kind the sum of the geometric phases vanishes. 
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